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High amplification events (HAEs) are common phenomena in ex- 
tragalactic gravitational lens systems (GLSs), where the multi- 
ple images of a distant quasar are observed through a foreground 
galaxy. There is a considerable brightness magnification in one 
of the quasar images during HAE. Grieger, Kayser, and Refs- 
dal (1988) proposed to use HAEs to study the central regions 
of quasars in GLSs. In this paper, we consider some problems 
concerning the identification of different source types on the basis 
of the HAE observations. We compare the results of light curve 
simulations to estimate a feasibility to distinguish different source 
models in GLSs. Analytic approximation methods yielding so- 
lutions of the lens equation in a vicinity of fold caustic crossing 
events are presented. The results are used to obtain amplifica- 
tion factors, which the higher-order corrections for the Gaussian, 
power-law, and limb-darkening models of a source take into ac- 
count. 



1. Introduction 

Gravitational lens systems can be viewed as natural tele- 
scopes that provide a valuable information about remote 
objects. In an extragalactic gravitational lens system 
(GLS), a distant quasar is observed through a foreground 
galaxy. The gravitational field of the galaxy can bend 
light rays from the source sufficiently so that there are 
multiple light rays that reach an observer. The observer 
sees an image in the direction of each ray, so that the 
source appears multiply imaged. 

Light rays from a quasar pass through the lensing 
galaxy in different regions. Local variations of gravi- 
tational fields in these regions, which are mainly due to 
a relative motion of the lensing galaxy and the source, 
lead to considerable brightness variations in each image 
(gravitational microlensing) , which can be detected even 



by modest telescopes. The typical time scales of these 
microlensing processes may vary from weeks to months. 
A comparison of the independent brightness variations 
in different images provides a valuable information about 
the lensing galaxy and about the source structure 
One of the most important applications of the microlens- 
ing concerns a unique possibility to study a fine structure 
of the central quasar region that cannot be resolved in 
another way with modern observational techniques. This 
is important because the quasars in GLSs have high red- 
shifts. Therefore, when studying the quasars, we learn 
something about a corresponding early epoch. 

HAEs - considerable brightness magnifications in one 
of the images of a quasar - are common phenomena in 
extragalactic GLSs. Grieger, Kayser and Refsdal [3] pro- 
posed to use HAEs to study the central regions of the 
quasars in GLSs. Typically corresponding variations of 
the brightness in a neighborhood of the HAE can be 
approximately described by a formula containing a few 
fitting parameters. This makes it possible to estimate 
some GLS parameters, in particular, the source size 0|. 
For example, in case of the well-known GLS Q2237+0305 
"Einstein Cross" several HAEs were observed [B^j 
and the estimates of the source size have been obtained 
for different source models [i-14|. Almost all HAEs in 
the Q2237+0305 GLS are attributed to the intersection 
of a fold caustic in the source plane (see, e.g., [Hi). In 
view of the increasing accuracy of photometric observa- 
tions, a possibility to distinguish different source models 
is also discussed (see, e.g., |16.[l7|). 

Here, we study some problems dealing mainly with the 
investigation of the source structure in an extragalactic 
GLS using the light curves of the source images. We dis- 
cuss the results of simulations of light curves for different 
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microlensed source models and present some analytic ap- 
proximations yielding solutions of the lens equation in a 
vicinity of the fold caustic crossing events. The results 
are used to obtain amplification factors for some source 
models that take the "post-linear" corrections into ac- 
count. 

Note that, in reality, we have, indeed, a unique light 
curve for each image of the quasar in a GLS. The high 
amplification events are not too frequent, and it may 
take a considerable time to accumulate a sufficient statis- 
tics and even to wait a repetition of the HAE. Thus, the 
main attention is paid to the determination of source 
properties from a single light curve compared with the 
light curves of the other images. On the other hand, 
since the work by Kochanek [18] followed by a number of 



works |15l . ll7l . ll9l - l23| , numerous statistical methods have 
been developed to process complete light curves. These 
methods treat the whole available light curves from all 
the images, not only HAEs. There is an enormous num- 
ber of the degrees of freedom which prohibit the simul- 
taneous determination of all the parameters of a GLS 
(microlens positions, their masses, and the source pa- 
rameters). But the observational data restrict possible 
realizations of the microlens masses and the positions, 
so one may estimate the conditional probability of cer- 
tain parameter values. Such approach is very attrac- 
tive because it allows one to use the whole aggregate 
of observational data on image light curves. However, 
this involves a large number of realizations of the mi- 
crolensing field and, respectively, a considerable com- 
puter time. 

On the other hand, the source structure manifests it- 
self only in HAEs; far from the caustic, the source looks 
like a point one, and all the information about its struc- 
ture is lost. If we restrict ourselves to a HAE neighbor- 
hood, then we use the most general model concerning a 
microlensing field described by a small collection of the 
Taylor expansion coefficients in the lens mapping. 

The structure of this paper is as follows. Section [2] 
contains some basics of the gravitational lensing, a short 
review of the problems concerning the determination of 
a source model from HAE observations, and a list of 
typical models. In Section [3l we present the results of 
simulations for these source models, in particular, the 
differences of light curves. In Section 21 we outline ap- 
proximation methods to investigate solutions of the lens 
equation used to derive the amplification coefficient of 
a point source. The result is used to obtain amplifica- 
tions of an extended source near the fold caustic for some 
source models. Section |6] contains the discussion of the 
results. 



2. General Relations and Notations 
2.1. The lens equation 

First, we recall some general notions concerning the 
gravitational lensing that can be found, e.g., in book 
[1]. The normalized lens equation has the form 



V$(x) 



(1) 



where <I> (x) is the lens potential. This equation relates 
every point x {xi, X2) of the image plane to the point 
y = {viiV2) of the source plane. In the general case, 
there are several solutions X(i) (y) of the lens equation 
([Ij that represent images of one point source at y; we 
denote the solution number by the index in parentheses. 

If there is no continuous matter on the line of sight, 
the potential must be a harmonic function: A<I> = 0. 
Below, we assume that this condition is fulfilled in a 
neighborhood of the critical point. We note, however, 
that if the continuous matter density is supposedly con- 
stant during a HAE, this can be taken into account by 
a suitable renormalization of the variables. 

The amplification of a separate image of the point 



source is 



(y) = l/|j(X(,) (y))| 



(2) 



where J (x) = \D {y)/D (x)| is the Jacobian of the lens 
mapping. In the microlensing processes, microimages 
cannot be observed separately; therefore, we need the 
total amplification that is a sum of the amplification 
factors of all the images. 

The critical curves of the lens mapping ([T]) are de- 
termined by the equation J (x) = and are mapped 
onto the caustic in the source plane. The stable critical 
points of a two-dimensional mapping can be folds and 
cusps only, the folds being more probable in a HAE. In 
this paper, we confine ourselves to the consideration of a 
fold caustic. When the point source approaches the fold 
caustic from its convex side, two of its images approach 
the critical curve, and their amplification tends to infin- 
ity. They disappear when the source crosses the caustic. 
These two images are called critical. 

2.2. Problems 

Let /(y) be the initial surface brightness distribution of 
an extended source. If the source center is located at the 
point Y = (11,1^2) in the source plane, then the total 
microlensed flux from the source is 



F(Y) 



/(y(x) - Y) dxidx2, 



(3) 
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X = {xi,X2)- The result of using Eq. ^ is obviously 
equivalent to the result of the well-known ray-tracing 
method f]\ (when the pixel sizes tend to zero). 
An equivalent representation of this formula is 



F(Y)= // K{y)I{y-Y)dyidy2, 



(4) 



where the point source amplification K{y) = "^Kiiy) 

i 

is the sum of amplifications of all the images. 

Near the caustic, one can approximate K{y) = Kq + 
Kcr{y), where Kq is the amplification of all noncritical 
images that is supposed to be constant during HAE, and 
Kcr is the amplification of the critical images. Due to 
a relative motion of the lensing galaxy and the source 
(quasar), the flux is a function of time representing the 
light curve of some quasar image in the GLS. In the 
lensing galaxy rest frame, the quasar motion may be 
considered as a straight-line motion with a sufficient ac- 
curacy. Here, we are mostly interested in the caustic 
crossing events when the quasar intersects the caustic in 
the source plane leading to a considerable enhancement 
of the image brightness. The light curve yields informa- 
tion about the source structure and about the gravita- 
tional field of the lensing galaxy. Particularly, the form 
of a light curve near a HAE depends on the function 
/(y), and the question is whether it is possible to use 
this dependence in order to identify this function. 

To illustrate the typical problems arising in the con- 
ventional treatment of a HAE, we consider an approx- 
imate relation for the amplification coefficient Kc^iy) 
of a point source located at the point y = {y 1,1/2) of 
the source plane which is located near the fold. In 
the appropriate coordinate system, the amplification of 
a point source can be approximated as i^^cr(yi, 2/2) = 



-1/2 



9(2/2) [1], where 2/2 is the distance between 



the point source and the caustic, and 6(2/2) is the Heav- 
iside step function. The observed radiation flux from the 
extended source image during the HAE is then obtained 
from Eq.(|4]) yielding an integral equation for the one- 
dimensional luminosity profile /(2/2) = / ^(2/17 2/2)^2/1' 

F{t) ^Ci+Ao f e(2/2)(2/2)-'/V(2/2 - Y2it))dy2, (5) 



Y(t) — {Yi{t),Y2{t)) is the source center trajectory 
which can be written as a linear function of time; F{t) is 
known from observations, Ci describes a contribution of 
noncritical images and can be considered to be constant 
during the HAE. Thus, we have an equation for f{y). 

The main problems concerning with this equation are 
as follows. 



(i) Equation ([5|) gives us nothing about the whole func- 
tion I{yi,y2) we are interested in, unless some supposi- 
tions about the form of the source are made, e.g., the 
supposition about circular symmetry. Also, even for 
known F{t) and Y2{t), Eq. ([5|) presents a kind of the 
ill-posed mathematical problems: small variations of in- 
put data can lead to a considerable change of the so- 
lution. A standard way to relax this difficulty involves 
additional restrictions on f{y) and/or the use of some 
explicit models for the brightness distribution /(y) con- 
taining a small number of free parameters. Some of these 
models are considered below in the next subsection. 

(ii) It is clear that a real brightness profile of the cen- 
tral quasar region is quite different from the simplified 
brightness distributions of the following section that can 
be considered rather as some reference models. How- 
ever, in view of the present-day accuracy of observa- 
tions, sometimes it can be difficult to distinguish even 
these simple source models using a HAE. For example, 
the authors of jl3| argue that the light curves from a ac- 
cretion disk can be well fitted with any brightness profile 
(Gaussian, uniform, etc.) of an appropriate source size. 
On the other hand, a number of authors (S- lol. 15- lol 23| 
discussed the delicate questions concerning the determi- 
nation a fine quasar structure from HAEs. For example, 
the authors of [l^l wrote that the GLITP data Q on 
Q2237+0305 admit only accretion disc models (see also 
iSlll)- Obviously, the presence of an accretion disk in 



the central region of a quasar is beyond any doubts, as 
well as the fact that the real quasar core can be quite 
different from its simplified models in question. The 
question is how to prove a posteriori the existence of 
the accretion disk on a basis of available observational 
data. 

One must also have in mind that, without using an 
additional information (besides HAEs), we cannot de- 
termine even the source size, because we do not know 
the value and the direction of the source velocity with 
respect to the caustic. Also, different ellipticities and 
orientations of the source may lead to different forms of 
the light curves during a HAE. 

(iii) The kernel Kcr of the integral equation ([SJ is a 
result of the so-called linear caustic approximation of the 
lens equation. This will work correctly only in the case 
where the source size is much smaller than the caustic 
curvature radius. Below, we consider corrections to 
that arise in the post-linear approximation. However, 
these approximations also require the source to be suffi- 
ciently small. This requirement can be violated in case 
of a complicated caustic network or in the presence of a 
population of small microlenses (cf. planetary masses) . 
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2.3. The extended source models 

Below, we list the simplest most commonly used bright- 
ness distributions of the source in a GLS; without loss 
of generality, they are chosen to be normalized to 1: 



JJ Iiy)dyidy2 = 1. 



To compare different models of brightness distribution, 
we have to use the same parameter that characterizes 
the size of an object. The most general is the r.m.s. size 

Rims • 



(6) 



For a slowly decreasing brightness profile (e.g., /(y) ^ 
|y| ck < 4 ) the r.m.s. size diverges. In the case of the 
circularly symmetric sources, the half-brightness radius 
Ri/2 is also widely used; it is defined by the relation 



Rl/2 

In J I{r)rdr 





oo 

vr J I{r)r dr. 




In the case of the Gaussian source model, 
1 



lair) 



exp 



f- 



(7) 



(8) 



where R — i?rins, R1/2 = i?-\/ln(2), and R stands for a 
size parameter. 

The limb-darkening (LD) model (see, e.g., yields 



where 



7ri?2 



E{r/R;q), 



(9) 



2\q 



s(C;9) = e(i-a(i-e 



and i?2 = (g + 2) i?2^jj,. Here, we assume q > 0. The 
half-brightness radius is R1/2 = Ry^l - (l/2)i/(9+i). 

Models ^ and ^ describe a class of compact sources 
with a fast decrease in the brightness. On the contrary, 
the power-law models describe a slower decrease 

at large r: 



7ri?2 



[1 + rVi?']' 



(10) 



where p > 1 is the exponent, and R is related to the 
r.m.s. radius i?rms as i?^ = (p — 2)i?2^g. Model (fTOj) 
is an alternative to The half-brightness radius of 



the source within this model is R1/2 = i?\/2i/(p-i) - l. 
For fixed i?rms and p — > 00, the brightness distribution 
tends to the Gaussian one. For small p, we have a 
"long-range" distribution; i?rms diverges for p < 2. 

Linear combinations of different distributions with dif- 
ferent parameters yield rather a wide class of symmetric 
source models. 

A more complicated profile is presented by the 
Shakura-Sunyaev model (AD) of accretion disk [25| 



lAuir) 



me{r - R) 



1 - 



(11) 



For the AD model, the half-brightness radius is R1/2 = 
AR, i?rms = 00 

Another model (l8l . [25| (ADl) can be defined as 



/adi(?') 



C 



API 
i?2 



exp(p3/4) _ 1 



p = r/R, (12) 



Cadi - 3(8^r(8/3)C(8/3))-i « 0.06, 

C{x) is the Riemann zeta function, T{x) is the Gamma 
function, and 



02 ^ „2r(16/3)C(16/3) 



r(8/3)C(8/3) 



21.4i?2 



It is also interesting to study a superposition of the LD 
and AD models (LA), where the accretion disk has a 
boundary 



/LA(r) = 



i?2 [exp(p3/4) _ 1] 1 



P = r/R, 



(13) 



where C(g) is a normalizing coefficient; C(l) = 
0.42, C(2) = 0.54, C(3) = 0.64. 

3. Simulations of Light Curves 

3.1. Equations of microlensing 

In the case without continuous matter on the line of 
sight, the lens equation takes the form 



N 



X - X,; 



(14) 



where are positions of the microlenses in the lens 
plane, and Re is the radius of the Einstein ring for one 
microlens which is assumed to be the same. Re = 1, for 
all microlenses. 
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1,0 -^ 




-3-2-10123 ' 



Fig. 1. Magnification pattern with the trajectory of the source 
in the source plane (coordinates are measured in units). The 
source moves from left to right along the straight line with uniform 
speed 

Here, we present the results of straightforward calcula- 
tions of the microlensed flux ^ to obtain the light curves 
for different realizations of the microlens positions. All 
calculations were performed for the microlensing optical 
depth <7 = 0.3. The total number of microlenses was 
1470. 

The microlens positions were chosen in a random way 
with uniform distribution over the field. The trajectory 
length has been taken so as to provide the caustic cross- 
ings. The size of the microlens field was chosen large 
enough to avoid boundary effects. 

3.2. Light curves 

The form of the light curve of the microlensed source 
in caustic crossing events depends on details of the 
source structure. However, the point is to determine 
the most relevant source model with regard for this 
form. In this work, we compare the light curves 
within the Gaussian, PL, LD, AD, and ADl source 
models. The simulations were performed for the set 
of 100 realizations of a microlensing field with the 

Table 1. Parameters of simulation 



Parameter 



Value 



Number of pixels 
Pixel size 

Source trajectory length 
Radius of field 

Microlensing optical depth (a) 
Source speed {V) 
Time discretization (St) 



1.23 X 10'^ 
0.01 Re 

2 -/?Ei 5 

70 Re 
0.3 
1 

0.01 Re/V 



Fig. 2. "Light curves": the amplification as a function of time 
for different source models that correspond to the magnification 
pattern in Fig. [T] 



optical depth a = 
of Q2237+0305 [I 



0.3 corresponding to parameters 
29|. All models have the same 



half-brightness radius Ri/2- The calculations were per- 
formed for the same microlenses fields. The source 
speed is = 1, so we can identify the source posi- 
tion as a function of time t. All the microlenses are 
static. The typical magnification pattern is shown in 

Fig.m 

First, we present the results of simulations with the 
same half-brightness radius R1/2 = 0.21; the power-law 
exponent was p = 3/2 for the "long range" PL model; the 
AD model also corresponds to this class of the power- 
law asymptotic dependence with p = 3/2). For LD and 
LA models, we have chosen q — 1 throughout the paper. 
From the "light curves" in Fig. [51 we observe a significant 
difference between the compact (LD and Gaussian) and 
"long-range" models. The long-range character of the 
latter reveals itself even at considerable distances from 
the caustic, where we expect that the brightness of all 
sources must have the same behavior as that of a point 
source. The differences between these two groups of 
models are essentially larger than the differences within 
each group (e.g., between the Gaussian and LD models). 
This conclusion is confirmed by the results of statistical 
considerations over 100 realizations shown in Table [5] 
for the half-brightness radius R1/2 — 0.21 as an exam- 
ple. 

To compare different models i and j, we used the rel- 
ative difference 



,^2maxl'^f;-^^fi' 



(15) 
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1,0 1,5 2,0 2,5 




Fig. 3. Magnification pattern with simple caustic crossing events 

where Ki{t) is the amphficatioii for the j-th model along 
the trajectory of source's linear movement. 

These results can depend on a complexity of the caus- 
tic involved into our consideration (i.e., there can be 
parts of the fold caustic close to the cusp points, or 
there can be dense aggregations of caustics). From 
many- year observations of the light curves of such GLS 
as Q2237+030, one can rule out such complex cases. 
Therefore, we considered some modification of our sta- 
tistical consideration with rather simple fold caustic 
crossings. However, the numerical results of this mod- 
ification with simple caustic crossing events (such as 
in Fig. I3|4I as an example) appeared to be nearly the 
same as those of Table [2j As an example. Table [3] 
shows the results of simulations analogous to those in 
Table [5] with the r.m.s. radius i?ims = 0.21; here, 
we excluded the models with i?ims = oo (here, p — 
3 for the PL model). The larger error is due to a 
smaller number of realizations with the "simple" caus- 
tic. 

Table 2. Relative difference between models in HAE 



i-th model 


j'-th model 


•n 


AD 


Gaussian 


0.074 ± 0.0012 


AD 


ADl 


0.085 ± 0.002 


AD 


LD 


0.091 ± 0.002 


AD 


PL 


0.038 ±0.001 


Gaussian 


ADl 


0.073 ±0.0017 


Gaussian 


LD 


0.042 ±0.001 


Gaussian 


PL 


0.073 ± 0.0013 


ADl 


LD 


0.094 ±0.002 


ADl 


PL 


0.052 ±0.0012 


LD 


PL 


0.090 ± 0.002 



Fig. 4. Light curves (magnification curves) for different models 
of source brightness profile that correspond to the magnification 
pattern of Fig. [3] 

3.3. Gaussian fittings of the accretion disk and 
limb darkening models 

The above results concern with a comparison of differ- 
ent models with the same R1/2 or i?rms- However, in 
reality, we do not know the source model to fit, and one 
may ask why a light curve is not fitted within a differ- 
ent model. Therefore, one must verify whether we can 
replace one model with a different one with some other 
source parameters to get a better fitting. 

We have fitted the limb-darkening and accretion disk 
model light curves with that of the Gaussian source of 
different radii. The half-brightness radius varied from 
i?i/2 = 0.2Re to i?i/2 = 0.24i?E; here, p = 3/2 and 
q = l. 

As we can see from Fig. [S] and Table 21 the Gaus- 
sian source cannot reproduce all the models, though the 
fitting results are rather good for the class of compact 
models. For example, the LD model can be replaced by 
the Gaussian model with the other source size. The rel- 
ative differences between models in Fig. [5] and Table 2] 



Table 3. Differences between source's models with the 
same iirms = 0.21R^ parameter in simple caustic crossing 
events; here, p = 3, q = 1 



i-th model 


j-th model 


V 


ADl 


Gaussian 


0.12 ± 0.04 


LD 


Gaussian 


0.03 ± 0.01 


PL 


Gaussian 


0.05 ±0.012 


ADl 


LD 


0.14 ±0.05 


ADl 


Power-law 


0.08 ±0.03 


LD 


PL 


0.08 ±0.02 
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-LA 
-LD 
AD1 

- Gaussian (R,,^^ 0.235 R^) 

- Gaussian (R,,j,= 0.23 R^ ... 0.205 RJ 
-Gaussian {R^^ = 0.2 R^^) 



1,0 



- LD-Gaussian (R„j=0.239 R^^) 
■ AD1-Gaussian (R„j=0.222 R,,) 
LA-Gaussian (R,,2=0.21 1 R^) 




Fig. 5. Light curves of the LD, ADl and AL models are fitted 
with a Gaussian source of different sizes from R1/2 = 0.2i?E to 
Ri/2 = 0.24i?E with step SR1/2 = 0.001i?E- These fight curves 
correspond to the magnification pattern of Fig. |3] 

show that it is possible to fit a considerable part of the 
whole curve with different source models on the accuracy 
level which is comparable to that of modern photometric 
observations. 



4. Approximation Methods 

Statement of the problem 

The study of a caustic crossing event is closely related 
to the investigation of the lens equation solutions. This 
equation near a fold can be expanded in powers of lo- 
cal coordinates; in the lowest orders of this expansion, 
the caustic is represented by a straight line; so, this ap- 
proximation is often referred as a "linear caustic approx- 
imation". In this approximation, the point source flux 
amplification is given by a simple formula ^ including 
the distance to the caustic and two fitting parameters. 
In the most known cases, the linear caustic approxima- 
tion is sufficient to fit the observed light curves during 
HAEs at the modern accuracy of flux measurements. At 
the same time, the consideration of "post-linear" terms 
is sometimes necessary to explain the present observa- 
tional data. The need for a modification of this formula 
- e.g., by taking the caustic curvature into account - is 
being discussed for a long time j^,[33,|3H. The correc- 

Table 4. Differences between models for fitted curves 



i-th model 


j-th model 


V 


AD 


Gaussian (/?i/2 = 


0.2) 


0.07 ±0.006 


LD 


Gaussian (R1/2 = 


0.24) 


0.026 ± 0.002 


LA 


Gaussian {R1/2 = 


0.2) 


0.034 ± 0.006 



Fig. 6. Relative differences between the Gaussian, LD, ADl, and 
LA models for best fitted curves 

tions to the amplification coefficient in the case of the 
macrolensing were the subject of investigations dealing 
with the problem of "anomalous flux ratios" [33 ■ A mod- 
ification of the post-linear approximation allows one to 
improve the quality of the fitting of the HAE light curve 
of the image C in the GLS "Einstein Cross" [33,34]. One 
may hope for that an improvement of the photometric 
accuracy of GLS observations will make it possible to 
obtain additional parameters of the lens mapping con- 
nected with the mass distribution in the lensing galaxy. 
Below, we study some points of an approximate solution 
of the lens equation near the fold caustic which probably 
is crossed by the source in many HAEs. 

The standard consideration of the caustic crossing 
events deals with the Taylor expansion of the potential 
near some point of the critical curve in the image plane. 
Let this point be the coordinate origin in an appropri- 
ate coordinate system, and let this point map onto the 
coordinate origin of the source plane. Further, we ro- 
tate the coordinate systems in both planes so that the 
abscissa axis on the source plane is tangent to the caus- 
tic at the origin, the quantity I2/2I defines locally the 
distance to the caustic, and yi defines a displacement 
along the tangent. For the harmonic potential, we can 
write the corresponding lens equation as 



yi = 2xi + a — X2) + 26x1X2 + c (xf — Sxix^) 
-d - 3x2x1) + gx\ + 

7/2=6 {x\ — X2) — 2axiX2 + d {x\ — 3x1X2) -l- 
-fc (xf-3x2X^) + fx\ + 



(16) 
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where a, b, c, d, g, and / are expansion coefficients. If the 
2/2 axis is directed toward the convexity of the caustic, 
then 6 < (at fold points, b ^ 0). 

We now proceed to the derivation of approximate so- 
lutions of Eqs. (IT6l) . To do this, we present two differ- 
ent methods |33|]. The first method deals with analytic 
expansions in powers of a small parameter. However, it 
results in nonanalytic functions of coordinates leading to 
nonintegrable terms in the amplification coefficient. The 
second method does not lead to such problems, though it 
uses a somewhat more complicated representation of the 
solution of the lens equation containing square roots of 
analytic functions. The methods agree with each other 
in a common domain of validity; moreover, we use the 
second method to justify some expressions in the ampli- 
fication formulas in terms of distributions to validate the 
applications to extended source models. A more detailed 
presentation may be found in fsS*!. 

First, we suppose that the source and the caustic lie 
on different sides from the yi axis. Then, j/2 > 0, and 
we substitute 



yi = t'^yi, xi^t^xi, X2^tx2, 



(17) 



where i = 1,2, and t can be considered as a parameter 
of vicinity to the caustic. This is a formal substitution 
that makes easier the operations with different orders of 
the expansion. After calculations, we put t = \ and thus 
return to the initial variables y^. However, if we put yi to 
be constant with varying then this substitution allows 
us to study a local behavior of critical image trajectories; 
t = Q corresponds to crossing the caustic by a point 
source, when two critical images appear. As was shown 
in [s^, this allows one to look for solutions of Eqs. (fT6|) . 
by using the expansions of Xi in powers of t: 

Xl = Xio + Xiit + Xi2t^ + 



y2 — —bx2 + t{—2axiX2 + 0X2)+ 



-t^{bx1 — 3dxiX2 + fS;2)- 



(19) 



The substitution of expansions p8|) into ([T9|) allows us 
to determine all coefficients successively. For example, 
for the zero-order terms, we have 



^ {vi - ^2/2) , ^20 = sVm/ \b \ , 



(20) 



where e = ±1 determines two different critical solutions. 
This approximation yields the well-known formula for 
the amplification ([5]). 

The first-order approximation terms have been derived 
in [s^ilsSj. In microlensing observations, two critical im- 
ages cannot be resolved, so we need the total amplifi- 
cation coefficient of two critical images. However, the 
contributions of the order of ^ i appear to be cancelled 
in calculations of the total amplification. Therefore, to 
obtain a nontrivial correction to the zero-order amplifi- 
cation, higher order approximations should be involved. 

[illiii. 



Note 



These corrections have been derived in 
that the corresponding second-order terms contain ex- 
pressions nonanalytic in y2- Though the second-order 
corrections are expected to be small, they appear to be 
noticeable in some cases even in the analysis of the avail- 
able data on light curves in the Q2237+0305 GLS. 

4.2. Method 2 

The other approach to the construction of approximate 
solutions of the lens equation in a vicinity of the fold is 
described in [s^- This allows us to present the critical 
solutions of system ([T9| in the following form: 



Xl ^ p + ire\Jw 



X2 ^ X20 + i2lt + X22t'^ + ■■■■ (18) 

It should be stressed that the analyticity in t does not 
mean that the coefficients of expansions ([T8| will be an- 
alytic functions of coordinates jji in the source plane (see 
below) . 

In terms of the new variables p7)) , system ^TE\\ takes 
the form (up to the terms ^ t^) 

j/i = 2xi — ax\ + t{2bxiX2 — dx^^ 

I .2/ ~2 o ~ -2 I ~4n 

-|-r [axi — 6cx\X2 + gx,2)^ 



X2 — ts + e\/w, e = ±1. (21) 

After the substitution of Eqs. ((2T|) in the lens equation, 
the terms containin g in teger and half-integer powers of 
w can be separated [33]. After some algebra, this yields 
a system that is convenient for an iterative procedure re- 
sulting in analytic expansions for the functions p, r, s, w 
both in powers of t and yi at every step of iteration (33|. 
In js^l, such a solution was obtained up to the terms 
~ t^. If we expand ^yw in powers of t, then we immedi- 
ately have the solution in the form (fT8|). 
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4.3. Amplification of a point source 

The solutions of the lens equation are then used to de- 
rive the Jacobians of the lens mapping (for both images 
near the critical curve). According to ([2]), the value of 
yields the amplification of individual images. As 
we pointed out above, we need the total amplification of 
two critical images (the sum of two amplifications of the 
separate critical images). In the second-order approxi- 
mation (using the expansion up to the terms t^), this is 



K„ = 



1 e (2/2) 



1 + Py2 + Qyi 



Kyi 
4 2/2 



(22) 



where the constants P and Q are expressed via the Tay- 
lor expansion coefiicients from Eq.(IT]), and 



2I&I ' 

8(2/2) is the Heaviside step function. Note that k is the 
caustic curvature at the origin which enters explicitly 
into the amplification formula. The parameters P and 
Q are independent; the explicit formulae for them can be 
found in |33ll34| . However, this is not needed when we 
use Eq. ([22|) for fitting the observational data, because 
these constants are whatever considered as free fitting 
parameters. 

Formula yields an effective approximation for the 
point source amplification near the coordinate origin 
provided that 1/2 > 0, and j/2/yi is not too small (see 
the term containing k) . For a fixed source position, this 
can be satisfied always by an appropriate choice of the 
coordinate origin, so that the source will be situated al- 
most on a normal to the tangent to the caustic. 

If the source is on the caustic tangent or in the region 
between the caustic and the tangent, then formula (|22p 
does not give a good approximation to the point source 
amplification. Nevertheless, in the case of an extended 
source, we will show that result (|22l) can be used to ob- 
tain approximations to the amplification of this source 
even as it intersects the caustic. However, to do this, we 
need to redefine correctly the convolution of ([22]) with a 
brightness distribution. 

5. Amplification of Extended Sources 

5.1. Transition to extended sources 

We now present the results of studies of the point 
source amplification within extended source models. Let 
/(y) be a surface brightness distribution of an extended 



source. If the source center is located at the point 
Y = {Yi,Y2) in the source plane, then the total mi- 
crolensed fiux from the source is given by 

Formula ([22| contains the nonintegrable term ^ 
0(y2)(2/2)~^^^- Therefore, the question arises of how 
formula (|22|) can be used in the situation where the ex- 
tended source intersects a caustic and some part of the 
source is in the zone between the tangent and the caus- 
tic. In view of Section |4l it is evident that the men- 
tioned term is a result of the asymptotic expansion of 
the root \/y2 + Kyft^/2 + ... in the approximate solution 
PT|) . The direct usage of a solution in the form (PT|) for 
the calculation of the Jacobians of the lens mapping and 
then for the derivation of amplifications does not lead to 
any divergences, and any nonintegrable terms in Kcr do 
not arise. Nevertheless, it is convenient to have a repre- 
sentation of Kci in the form of an expansion in powers 
of a small parameter. Such an expansion can be carried 
out correctly after the substitution of Kct into integral 
(jl]). On this way, starting from ([2T|) . it is easy to show 
that, to define Kcr correctly, one must replace the term 
0(y2)(2/2)~^^^ in (|22|) by the distribution (generalized 



function) (1/2) 
-3/2 



-3/2 



36]. We recall that the distribution 



2/_|_ of the variable y is defined by the expression 



y+^^^fiy)dy 



-i/2^/(y) 

dy 



dy 



for any test function f{y). 
After this redefinition, we have 



0(2/2) 



2V\b\m 



[1 + Py2 + Qyi ] 



yl{y2) 



-3/2 



(23) 



This formula can be used to correctly derive an approx- 
imate amplification of a sufficiently smooth extended 
source including the case where the source crosses the 
caustic. 

5.2. Gaussian source 



Formula ([23|) has been used [33| to derive the amplifica- 
tion of a Gaussian source with the brightness distribu- 
tion ([5]), the limb-darkening source, and the power- law 
source (see the next subsections). 

Further, we use the dimensionless coordinates s = 
Yi/i?, h~Y2/Roi the source center and the functions 
00 

4 {h) = j u^-^l"^ exp {-u^ + 2u\i) du = 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 4 



397 



A.N. ALEXANDROV, V.M. SLIUSAR, V.I. ZHDANOV 



(24) 



71=0 



These functions can be expressed in terms of the conflu- 
ent hypergeometric function iFi or the parabohc cyhn- 
der function D: 



The substitution of and (|H]) in (g]) yields 



(25) 



R 



P$i (h) - (h) + Qs*o {h) - KsH2 {h)\ } 



Here, 

$0 {h) = Ioih)exp{-h^) , 
^i{h)^h{h)exp{-h^), 
$2 (h) = [hlo (h) - h m exp {-U") . 



Note that the main term of ([26|) which corresponds to 
the hnear caustic approximation was first obtained in 
work fSTj . 

5.3. Limb- darkening source 

Analogous considerations allow us to obtain formulas 
for the amplification of extended sources for the limb- 
darkening and power-law brightness profiles; the results 
are represented analytically in terms of the hypergeo- 
metric function 2-F1 [38|. 
Denote 



Xfc,, [h) = l^^^f, I V^--- S(y -h-q+ 1/2) dy, 



k=l,2. We have 



Xk,gih)^2''+Hi + hr- 



_,^, r(g + 2)r(fc + i) 

T{q + k + 2) 



X 2F1 



1 



k + 2: 



' 2'^ 2 
for — 1 < /i < 1 and 



for h> 1. 

For fc = — 1, we define 

X_i,g (h) = 4 (q + 1) (/lXo,,-i - Xi^,_i) . 

Then, in case of the model with limb darkening ([9]), 
the critical images disappear when the source lies on 
the outer side of the caustic (i.e., for h < —1). The 
amplification due to critical images takes the form 



(26) -^(^X_i.,+i(/i) + Q sXo,,{h) - J .2X_i.,(/i)] } 

5.4. Amplification for a power-law source 

The result for the amplification involves integrals: 



T{p-l) 



2 dy 



p-1/2 



= Y^"" + i 2p - fc - (1 + /.^)^/2+3/4-p ^ 

/, 1 ,31/ h 

''^^^['^r''-'-r'-'2['^VfW 

for k = 0,1, B {x, y) being the Beta-function. 

We extend this to fc = — 1 in view of the definition of 
iy)+^^'^- So, we have 

^-i^p{h) = 4(p - l)[h^o^p+i{h) - *i,p+i(/i)]. 

Now. the amplification due to critical images takes the 
form 

2\/7r &LR I- L 



-^-^*_i,p_i(/i) + Qs«'o.p(/i) - ^s2*_i 



pih)]} 



The zeroth approximation to this formula has been de- 
rived in 0. 
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6. Discussion 

Prior to sum up the results of this paper, we must 
stress that the consideration of HAEs without some ad- 
ditional information may only provide a very uncertain, 
on order-of-magnitude level, and even ambiguous infor- 
mation about the source parameters. For example, one 
may expect that, by adding a sufficient number of small 
microlenses, it is possible to reproduce fine features of 
a light curve during a HAE and to fit the light curve of 
any compact source. The introduction of planetary mass 
objects in the lensing galaxy (which is quite reasonable, 
see, e.g., ^39-42]) yields a fine caustic mesh, which, in 
turn, can give rise to fine features of light curves dur- 
ing a HAE similar to those due to any given brightness 
distribution over the source. Therefore, we point out 
the assumptions and the class of models involved, which 
is necessary for a reasonable formulation of the prob- 
lem. 

The simulations of the present paper use the GLS pa- 



We 



rameters similar to those of Q2237+030 
consider the equal mass microlensing system; we do not 
consider any mass distributions and/or populations of 
small (planetary) masses. Next, we consider the most 
simple source models without effects of ellipticity, etc. 
Nevertheless, we see that, even under these rather severe 
restrictions, there is an ambiguity in the determination 
of the source model from the observations of HAEs. The 
limb-darkening source can be successively fitted with the 
Gaussian source model (with the other R1/2) on 1-2% ac- 
curacy level. This confirms the conclusions of Morton- 
son et al. [17] (see also [43]) that the surface brightness 
profile has little effect on microlensing. What can be de- 
termined within the modern accuracy is the source size; 
one can also distinguish either we deal with a compact 
source or not (i.e., with a steep decrease of the brightness 
for large radii). 

The investigation of the source structure from HAEs 
is closely related to the solution of the lens equation in 
the caustic region. There are also different approaches 
to this problem under different restrictions. In this pa- 
per, we confined ourselves to the case of the fold caus- 
tic. We outlined two methods that enable us to obtain 
the critical solutions of the gravitational lens equation 
near a fold with any desired accuracy fss', "s^]. In or- 
der to obtain nontrivial corrections to Kd- obtained in 
the linear caustic approximation, the higher orders of 
the expansion of the lens equation must be taken into 
account, as compared to works [32|,|35|. The modified 
formula for contains three extra parameters. This 
is applied to the Gaussian, power-law, and limb dark- 



ening models of an extended source. The fitting of the 
light curve of GLS Q2237+0305C using these modified 
relations shows [33*, 'si'] that some of these corrections 
can be statistically significant even at the present accu- 
racy level. This means that if we are looking for some 
fine effects in HAEs due to the source size, a consistent 
treatment must involve sometimes the higher-order cor- 
rections to solutions of the lens equation. 

This work has been supported in part by the "Cos- 
momicrophysics" program of the National Academy of 
Sciences of Ukraine. We are grateful to the staff of 
the Computer Center of National Technical University 
of Ukraine "Kyiv Polytechnic Institute", where a consid- 
erable part of our simulations has been performed. 

1. P. Schneider, J. Ehlers, and E.E. Faiko, Gravitational 
Lenses (Springer, New York, 1992). 

2. J. Wambsganss, Gravitational Lensing: Strong, Weak, 
and Micro, edited by G. Meylan, P. North, and P. Jetzer 
(Springer, Berlin, 2006), p. 453. 

3. B. Grieger, R. Kayser, and S. Refsdal, A&A 194, 54 
(1988). 

4. J. Huchra, V. Gorenstein, S. Kent, I. Shapiro, G. Smith, 
E. Horine, and R. Perley, AJ 90, 691 (1985). 

5. D. Alcalde, E. Mediavilla, O. Moreau, J.A. Munoz, 
C. Libbrecht et al., ApJ 572, 729 (2002). 

6. A. Udalski et al. Acta Astron. 56, 293 (2006). 

7. P.R. Wozriiak, C. Alard, A. Udalski, M. Szymariski, 
M. Kubiak, G. Pietrzyriski, and K. Zebruri, ApJ 529, 
88 (2000). 

8. M.B. Bogdanov and A.M. Cherepashchuk, Astron. Repts. 
46, 626 (2002). 

9. V.N. Shalyapin, Astron. Lett. 27, 150 (2001). 

10. V.N. Shalyapin, L.J. Goicoechea, D. Alcalde, E. Medi- 
avilla, J.A. Munoz, and R. GiLMerino, ApJ 579, 127 
(2002). 

11. J.S. Wyithe, R.L. Webster, and E.L. Turner, MNRAS 
309, 261 (1999). 

12. J.S. Wyithe, R.L. Webster, and E.L. Turner, MNRAS 
318, 762 (2000). 

13. J.S. Wyithe, R.L. Webster, E.L. Turner, and D.J. Mort- 
lock, MNRAS 315, 62 (2000). 

14. A. Yonehara, A J 548, 127 (2001). 

15. R. Gil-Merino, J. Gonzalez-Cadelo, L.J. Goicoechea, 
V.N. Shalyapin, and G.F. Lewis, MNRAS 371, 1478 
(2006). 

16. L.J. Goicoechea, D. Alcalde, E. Mediavilla, and 
J.A. Muiioz, A& A 397, 517 (2003). 

17. M.J. Mortonson, P.L. Schechter, and J. Wambsganss, 
ApJ 628, 594 (2005). 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 4 



399 



A.N. ALEXANDROV, V.M. SLIUSAR, V.I. ZHDANOV 



18. C.S. Kochanek, ApJ 605, 58 (2004). 

19. T. Anguita, R.W. Schmidt, E.L. Turner, J. Wambsganss, 
R.L. Webster, K.A. Loomis, D. Long, and R. MacMillan, 
A&A 480, 327 (2008). 

20. S. Poindexter and C.S. Kochanek, ApJ 712, 658 (2010). 

21. S. Poindexter and C.S. Kochanek, ApJ 712, 668 (2010). 

22. S. Poindexter, N. Morgan, and C.S. Kochanek, ApJ 673, 
34 (2008). 

23. V.G. Vakulik, R.E. Schild, G.V. Smirnov, V.N. Dudinov, 
and V.S. Tsvetkova, MNRAS 382, 819 (2007). 

24. M. Dominik, MNRAS 353, 69 (2004). 

25. N.I. Shakura and R.A. Sunyaev, Astron. Astrophys. 24, 
337 (1973). 

26. R. Schmidt, R.L. Webster, and F.G. Lewis, MNRAS 295, 
488 (1998). 

27. X. Dai, E. Agol, M.W. Bautz, and CP. Garmire, ApJ 
589, 100 (2003). 

28. E.V. Fedorova, V.L Zhdanov, C. VignaU, and 
G.G.C. Palumbo, A&A 490, 989 (2008). 

29. V. Vakuhk, R. Schild, V. Dudinov, S. Nuritdinov, 
V. Tsvetkova, O. Burkhonov, and T. Akhunov, A&A 
447, 905 (2006). 

30. C.J. Fluke and R.L. Webster, MNRAS 302, 68 (1999). 

31. A.B. Congdon, C.R. Keeton, and C.E. Nordgren, MN- 
RAS 389, 398, (2008). 

32. C.R. Keeton, B.S. Gaudi, and A.O. Fetters, ApJ 635, 35 
(2005). 

33. A.N. Alexandrov and V.L Zhdanov, e-print 
larXiv: 1006.5903 1 (2010). 

34. A.N Alexandrov, V.I Zhdanov, and E.V. Fedorova, As- 
tron. Lett. 36, 329 (2010). 

35. A.N Alexandrov, V.I Zhdanov, and E.V. Fedorova, Visn. 
Kyiv. Univ., Astron. 39-40, 52 (2003). 

36. I.M. Gel'fand and G.E. Shilov, Generalized Functions, 
Vol. 1 (Academic Press, New York, 1964). 



37. P. Schneider and A. Weiss, A&A 171, 49 (1987). 

38. H. Bateman and A. Erdelyi, Higher Transcendental Func- 
tions, Vol. 1 (McGraw-Hill, New York, 1953). 

39. M. Jaroszynski, B. Paczynski, AA 52, 361 (2002). 

40. S. Mao and B. Paczynski, ApJ 374, 37 (1991). 

41. Y. Tsapras, K. Home, S. Kane, and K. Carson, MNRAS 
343, 1131 (2003). 

42. J.S.B. Wyithe and E.L. Turner, MNRAS 320, 21 (2001). 

43. A.B. Congdon, C.R. Keeton, and S.J. Osmer, MNRAS 
376, 263 (2007). 

Received 30.01.11 

EOEKTH HEPETHHy KAYCTHKH I MO/JEJII /],>KEPEJIA 
y rPABITAIJIHHO-JlIHSOBHX CHCTEMAX 

O.M. AjieKcandpoe , B.M. Cnmcap, B.I. }KdaHoe 
P e 3 K) M e 

no.z],ii 3 BejiHKHM ni^CHjieHHHM (nBn) e 3BHMaHHHM HBHmeM y no- 
3arajiaKTHHHiix rpaBiTari;iHH0-jiiH30BHx CHCxeMax (FJIC), fle cno- 
CTepiraiOTb ^eKijiBKa 3o6pa?KeHB Bi^ajienoro KBasapa Ha 4)OHi ra- 
jiaKTHKH, mo 3Haxo;i,HTbCH Ha nepe^ii^HBOMy Hjiani. UpoTHroM riBII 
Bifl5yBaeTbca SHanne 36ijibmeHHa acKpaBOCTi b oflHOMy la 3o5pa- 
JKCHB KBa3apa. Fpirep, KaH3ep Ta Pe4)Cflaji 3anponoHyBajiH bhko- 
piiCTOByBaTH nBn AJia BiiB-^eHHa i^eHTpajiBHHx oGjiacTefi KBa3apa 
B rjlC. y uiii CTaTTi ^ocjii^>KyeMO noB'H3aHi 3 u,hm HHTanHa, mo 
CTOcyiOTBca ifleHTH4)iKau,ii pi3HHx thhIb fl>Kepejia na 5a3i cnocTe- 
pe^KCHB nBn. Mil HopiBHioeMO pe3yjii>TaTn hhcjiobhx MO.zi,ejiiOBaHB 

KpHBHX 6jIHCKy A^ia TOrO, III,o6 OII,iHHTH MO>KJlHBiCTb Biflpi3HHTH 

MO.zi,ejii .njKepejia b FJIC. SanponoHOBano Ha6jiH>KeHi cxeMH AJia 
po3B'a3aHHa piBHanna jiiH3H b okojiI KayCTHKH — CKJia^KH. Pe- 
3yjibTaTn BHKopHCTaHO ^ji^ OTpHMaHHa Koe<|)iii,ieHTiB niflCHjieHHa, 
mo BpaxoByioTb nonpaBKH bhcokofo nopaflKy fljia MOflejieii rayclB- 
CbKoro i CTeneneBoro fljKepejia, a TaKOjK fljKepejia 3 noTeMHiHHaM 
flo Kpaio. 
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